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Abstract 

We suggest a simple synthesis model of an eclipsing binary system 
which includes one component with strong stellar wind. Numerical sim- 
ulations show that the shape of the light curve (and in particularly the 
widths of the minima) strongly depends on wind parameters. Wind effects 
are crucial in modelling light curves of binaries including e.g., WR stars. 

1 Introduction 

Several very massive Wolf-Rayet (WR) stars were recently discovered in 
eclipsing binaries, e.g. WR20a with WR mass Mwr ~ 80 Mq (Bonanas et 
al., 2004), NGC3601-A1 with Mwai « 116 M©, Mwr2 « 89 Mq (Schnurr 
et al., 2008), etc. The radii of WR stars in these systems turned out to be 
surprisingly large (~ 2QRq) (Bonanas et al., 2004; Moffat et al., 2004), 
similar to radii of main sequence (MS) stars of comparable masses. Ac- 
cording to common evolutionary scenarios, a WR star is formed through 
a massive loss of its hydrogen-rich outer layers. One could expect that 
resulting stellar radii should be smaller than radii of MS stars with com- 
parable masses. 

The radii of the WR stars in the above binaries were obtained by 
analysing eclipses. Absorption in extended winds of WR stars was not 
accouned for, instead a standard Roche model (Wilson, 1979) of stars with 
thin atmospheres was used for both components of a binary. Including 
extended wind in a model is expected to reduce the radii and hopefully 
make them compatible with predictions of theory. Earlier attempts to 
account for wind absorption were made in Pustyl'nik & Einasto (1984), 
Antokhina & Cherepashchuk (1988), Kallrath & Milone (1999). In these 
papers, it was assumed that wind expands with constant velocity. The 
only advantage of such an assumption is simplicity of calculations. At the 
same time it is clearly inadequate to actual situation in hot star winds 
expanding with acceleration. In current paper, we present a simple wind 
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model which allows one to use any parametric form of the velocity law for 
the wind. 



2 Model 



Roche model (Antokhina, 1988; Antokhina et al., 2000) is used for both 

components. One component is surrounded by homogeneous, spherically 
symmetric wind. While computing the model light curve, emission from 
every elementary area of each Roche surface is summed up, accounting for 
eclipses by stellar bodies and for absorption in the wind. Velocity law is 
set by a parametric expression, e.g. the /3-law: v{r) = Vao (l — ^r)''- The 
main opacity in optical continuum is electron scattering, optical depth 
from the point {p,zo) in the wind (see Fig.l) is 



t{p,zo) = J e{z)dz 

^0 



(1) 



where the opacity coefficient 



e(r) = (TTne('r) = 



(2) 



47rmpv(r)r2 ' 

Ze is the number of electrons per unit atomic mass. If the chemical abun- 
dance of the wind is solar, hydrogen being the most abundant element, 
«e = 1 (in hot star winds hydrogen is ionized). In WR winds the most 
abundant element is helium, so Ze ^ 0.5 if helium is fully ionized, Ze ^ 0.25 
in the Hell zone. 

The model parameters along with their description are listed in Ta- 
ble 1. 




Observer 



Figure 1: The coordinate system for the wind model 
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Table 1. Input parameters of the model. 



Parameters Description 



Binary parameters 

Mi,M2 (Mq) Stellar masses 

P (d) Period 

e Eccentricity 

Lo Longitude of periastron, star No.l 

i (deg) Orbital inclination 

m, iJ,2 Roche lobe filling coefficients (/i = 1 for star filling 
its inner critical Roche lobe at periastron position) 

Ti,T2 (K) Average effective temperatures of the components 

/3i ,02 Gravity darkening coefficients 

A\,A2 Bolometric albedos 

F\ , F2 Ratio of surface rotation rate to synchronous rate 

x\, X2 Limb darkening coefficients 

Z3 Third light 

A(n) Effective wavelengths of monochromatic light curves 

Wind parameters 

M (Mq I year) Mass loss rate 

Voo (km/sec) Terminal velocity of the wind 

Parameter of the velocity /3-law 

2e Number of electrons per unit atomic mass 



3 Results 

To illustrate the influence of wind opacity on model light curves, we com- 
puted two models, a detached binary (Model 1) and a contact binary 
(Model 2). In each model stellar parameters of both components are 
identical but the secondary component is surrounded by a wind. This 
configuration helps to isolate the effects of the wind opacity on the light 
curve. When the wind is absent, the two minima of the light curve are 
identical. When the wind is "turned on" , the depth of the secondary min- 
imum is dramatically increased, while the depth of the primary minimum 
is decreased (in the secondary minimum the component surrounded by the 
wind is in front). The stellar parameters in both models were fixed while 
wind parameters varied. The values of model parameters are listed in Ta- 
ble 2. Wind parameters were varied one by one. While varying one param- 
eter, the other two were fixed. The fixed values are M = 2-10"'^ Mq /year, 
Voo = 2000 km/s, /3 = 0.5. We assumed that wind consists of fully ionized 
helium {ze = 0.5). 

Simulations show that the shape of model light curves strongly depends 
on wind parameters (sec Figs. 3, 5). Clearly there is some degeneracy 
between the parameters. It is unlikely that the three most eccential wind 
parameters can be obtained independently from light curve analysis alone. 
Luckily, when solving light curves, most of them can be fixed at values 
obtained from independent sources (e.g., Voo from UV spectra, A4 from 
radio data etc.). In light curve analysis, the increase in the width of the 
minima most importantly influences radii and temperatures of stars. Our 
simulations show that it is absolutely crucial to take into account extended 
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Figure 2: Sky-plane view of Model 1 (left) and Model 2 (right) 



winds while modelling light curves of binary systems including WR stars. 
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Table 2. Parameters of simulated models. 



Parameter Model 1 (detached Model 2 (contact 

configuration) configuration) 



Binary parameters 

Ml (Mq) 10 10 

M2 (Mq) 10 10 

P (d) 5 5 

i (deg) 80 70 

til 0.5 1.0 

At2 0.5 1.0 

Ti{K) 30000 30000 

T2IK) 30000 30000 

Fi,F2 1.0 1.0 

/3i,/32 0.25 0.25 

Ai,A2 1.0 1.0 

XI, X2 1.0 1.0 

e 0. 0. 

LJ 0. 0. 

I3 0. 0. 

A (A) 4400 4400 

Wind parameters 

M (Mg/year) 10"^ - 4 ■ 10"^ 10"^ - 4 ■ 10"^ 

Voo (km/sec) 1000 - 2000 1000 - 2000 

/3 0.5-1.5 0.5-1.5 

0.5 0.5 
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Figure 3: Influence of mass loss rate M on the shape of light curve. Model 1 
(left), Model 2 (right). 




Figure 4: Same as in Fig. 3 , for the terminal wind velocity Vc 




Figure 5: Same as in Fig. 3, for the parameter /3 of the velocity law 
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